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We construct a unifying theory of geometric effects in mesoscopic stochastic kinetics. We demon- 
strate that the adiabatic pump and the reversible ratchet effects, as well as similar novel phenomena 
in other domains, such as in epidemiology, all follow from geometric phase contributions to the 
effective action in the stochastic path integral representation of the moment generating function of 
particle fluxes. The theory provides a universal technique for identification, prediction and calcula- 
tion of pump-like phenomena in an arbitrary mesoscopic stochastic framework. 



PACS numbers: 03.65.Vf,05.10.Gg,05.40.Ca 

Introduction. A number of effects in classical nonsta- 
tionary statistical physics, such as the reversible ratchet 
[U [21 13] and the adiabatic pump effects [41 15] , are known 
(or anticipated) to have a geometric origin. Examples 
and applications can be found in a variety of fields, in- 
cluding electronic metrology [B], cell motility [7j, ion 
pumping through a cell membrane [Hj, and dissipative 
chemical kinetics [9J . The distinct feature of these effects 
is that, under a slow periodic perturbation, transport co- 
efficients are not a simple average of those in the strict 
static approximation, but contain an extra component, 
which changes its sign under a time-reversal of the per- 
turbation. Although these effects have been well-studied 
in a variety of fields (see [10] for reviews) , a general the- 
ory, which clearly disambiguates the pump (or ratchet) 
currents from other nonequilibrium transport, provides 
a unified view of disparate pump-like phenomena, and 
suggests universal quantitative methods for calculation 
of moments of pump fluxes, is still missing. 

We address this problem using the recently intro- 
duced stochastic path integral representation of the mo- 
ment generating function (mgf) of fluxes in mesoscopic 
stochastic systems [T^l |T3]. In this letter, we demon- 
strate that the stochastic path integral technique can be 
employed to calculate moments of pump fluxes in a gen- 
eral stochastic driven system in the mesoscopic (many 
particles) and the adiabatic (slowly varying external pa- 
rameters) regimes, and that it makes a clear distinction 
between the pump fluxes and other currents by relating 
the former to a geometric phase contribution to the flux 
mgf. Our theory clarifies the connection between the re- 
versible ratchet and the adiabatic pump and allows to 
identify similar effects in new contexts, which we demon- 
strate in a specific scenario from epidemiology. 

Pump current from particle exclusion. Let two absorb- 
ing states S and P (substrate and product in a Michaelis- 
Menten enzymatic reaction [SJ [TT] , distinct cellular com- 
partments, neighborhoods of a city), exchange particles 
(molecules, humans) via an intermediate system B (bin, 
enzyme, membrane channel, transportation hub). Our 
goal is to find the S P flux J and its fluctuations 



on time scales much larger than the fluctuation time in 
the bin B, assuming the mesoscopic regime with a large 
typical number of particles in the bin, iV ^ 1. 

Particles interact, and the in- and out- going transition 
rates can depend on the number of particles in the bin, 
N{t). The simplest example of this kind is when the bin 
has a finite size, so that N < Nb — const < oo. Then the 
in-rates are proportional to the number of empty spaces 
in the bin, while the per-particle out-rates are not af- 
fected by the occupancy. The full kinetic scheme is 

1. S ^ B; rate ki{N,t) = qi{t){NB - N), 

2. B^ S; rate fc_i(iV) = q^^N, 

3. P^B- rate k^2{N,t) = q_2(<)(A^B - N), 

4. B ^ P; rate k2{N) = q2N. 

We allow <7i and g_2 to undergo a slow periodic modula- 
tion with a frequency to, which can be achieved in the bio- 
chemical context by coupling S and P to particle baths 
with modulated chemical potentials. In other transport 
problems, such as transportation systems, the same mod- 
ulation may be produced by time-of-day variations. We 
note that, unlike in jl2j , our formulation has three time 
scales: fast instantaneous jumps among states, equilibra- 
tion in the bin, and adiabatic changes of the rates. 

Now the path integral technique of [12] can be ap- 
plied. Since N ^ 1, there exists a time scale St, 
over which many transitions into and out of B hap- 
pen, but the fractional change in the bin occupancy re- 
mains small, 1 <C 5N <C N. Then the rate changes 
6ki, i = —2,-1,1,2 are also small, and all transitions 
are uncorrelated and Poissonian. Thus the probability 
of the number of particle transitions for the i'th re- 
action over time St, denoted by SQt, can be written 
as P{SQ,,t) = ±. J;^^dx^e-'^'^'^^+^^"^^^^''^^', where 
NsH^St = k^{N,t)[exp{ixi) - l]St = ki{N,t)e^.St is the 
mgf of a Poisson process with the mean kiSt. Note that 
we define = e*^ — 1 for any x. 

Our goal is to find the mgf of the net particle number 
Qp transfered into P during a long time interval (0,T). 
This is formally given by an integral over fiuxes at each 



2 



moment of (discretized) time weighted by probabilities 
P{5Qi,t) and constrained by particle conservation laws: 

{e'xcQp)^ Y{dN{tk) n d5Q,{tu)PmXtk)]x 

fc=l i=±2,±l 
X ^^Xc(SQ2(t,)~SQ-2(t^)) S[N{tk+i) - N{tk) 

- 5Qi{tk) - SQ-2{tk) + SQ^iih) + SQ2{tk)]. (1) 

Here we used the identity Qp = J2k=i[SQ2{tk) — 
SQ-2itk)], and we introduced a variable xci which is con- 
jugated to Qp and "counts" particle transfers into/out 
of P. Now, using the Fourier representation of the S- 
function, we integrate over 5Qi{t}S) and Xi(^fc) i^i Q ^^^id 
reduce it to a path integral over N and its conjugated x 



by qi{t). For the periodic driving, as we consider here, 
with a period Tq = iir/oj and with fixed g_i and q2, we 
rewrite the contour integral as the integral of Fq^^q_^ — 
dqi^-2^dq_^Ai over the surface Sc enclosed by c. Then 
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where 
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- (f A - dq^ ^ I dqidq^2Fq,,q_^{q), (9) 

Jc J^O JSc 



P,^,,-A^) m^q^e^^ + Aq_,q_2e-^a + K^?'^ 



(10) 



5ci = ^ dt{K- \lK^ + 4gig2exc + ^q-iq~2e 



xc ^ •*y-iy-2'i-xc 
(11) 



where all pre-factors are absorbed into the measure, and The 2-form ^^.^^^.(q) is an analog of the Berry curva- 
ture in quantum mechanics. As follows from ([9]), nonzero 
Berry curvature is responsible for the reversible compo- 
nent in the particle fluxes. Its presence in our model is 
due to particle exclusion within the central bin. If fci and 
fc_2 were independent of N, Fq^^q_^ would be zero. 

Now moments of the flux between absorbing states can 
be derived easily by differentiating ([s]) with respect to xc- 
In particular, the average flux is 



i7(x, iV, t) = [qi{t)e^^ + q.2{t)e-i^+^c)\ (1 " N/Nb) 
+ [(i-i^x + 92e(x+xc)] ^I^B- (3) 

The explicit time-dependence of H is due to the slow 
periodic changes in qi{t) and q-2{t)- 

The exponent in ^ has a factor of Nb in it. Thus for 
Nb — > oo, the path integral is dominated by the saddle 
point or classical values Xci and A^ch 



*^^' = — — ' 



dXci 



(4) 



J=Ju 



Jc] 



Here the boundary terms disappear, as explained in de- 
tail in [13]. Furthermore, since the Hamiltonian ([s]) is 
linear in N there are no higher order in 1/A'^ corrections. 

Since we assume adiabatic variation of qi and q_2- the 
quasi-equilibrium is a good approximation to the exact 
solution of Q . It corresponds to setting time-derivatives 
in Q to zero and treating time as a parameter. For 
periodic rates variations, this leads to 



92 + q-i 
dqidq-2 ^ ^3 
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2(gi + g-2e-*xc-) 
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qi + g-2e"*xc + (q_-^ + q^e'xc)e'2ixai 



(5) 
(6) 



where w denotes the accuracy of 0{uj/qi), K — qi + 

q-2 + q-i + q2, K- ^ qi + 9-2 - q-i - q2- Since the 
Hamiltonian is quadratic in its arguments near the saddle 
point, corrections of the order 0{uj/qi) in ([sj [6| lead to 
0[{uj/qi)'^] contributions to the mgf, setting the accuracy 
of our results. We now have 



where the pump term is due to the particle interactions 
and the corresponding geometric contribution, while the 
classical flux would exist even in the stationary limit. No- 
tice that this flux is Nb times its value for a single driven 
Michaelis-Menten enzyme, studied in [5]. The same holds 
for the entire mgf, and hence for all flux moments. Thus 
we refer the reader to [5] for further analysis of the model. 
Here we note that this scaling is not a coincidence since 
the current model is, indeed, equivalent to Nb indepen- 
dent enzymes, where N corresponds to the number of 
enzyme-substrate complexes. 

In [5] , we used an analogy with the quantum mechan- 
ical Berry phase to derive the pump flux (12 1, while a 



(7) 



where the vector A, Ai = iXci{dq-Nci) /Nb, is deflned in 
the space of parameters qi, and the contour c is given 



classical stochastic treatment, which does not require di- 
agonalizing large matrices, is being used now. Existence 
of these alternative approaches is not surprising because 
any discrete quantum mechanical system can be mapped 
onto a mathematically equivalent classical Hamiltonian 
system [14], and then the Berry phase transforms into 
a dynamic contribution to the classical action [TS]. The 
contribution of the present derivation is to show that one 
can derive the classical Hamiltonian for a discrete Markov 
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chain by considering many identical independent copies 
of the system and identifying the per-copy contribution 
by taking the large copy number limit. Alternatively, 
one can derive the classical reformulation directly from 
the Schrodinger's equation as well ^14. ,15. . 

The reversible ratchet effect. Now we show that the 
geometric contribution to the mgf is responsible also for 
the ratchet effect in a periodic potential. Consider a sys- 
tem of noninteracting particles moving in a periodic po- 
tential V{x,t), which changes adiabatically with time so 
that V{x,t) = V{x,t + To) and V{x,t) = V{x + L,t). 
In the overdamped case, the average density of particles 
satisfies the Fokker-Plank equation 

dtp{x,t) = -dx[A{x,t)p{x,t)] + Ddlpix,t), (13) 

where D is the diffusion coefficient, and A{x, t) = 
—dxV{x,t) is the force. The current in this model under 
an adiabatic deformation of the potential was previously 
studied in [5], and the similarity of the final expression 
and the Berry phase in quantum mechanics was pointed 
out. The close connection between the classical ratchets 
and the Berry phase also has been anticipated in [Ij [3] . 
In our following rederivation, we explicitly show that the 
ratchet current has its origins in the geometric phase. 
Namely, it emerges from the complex geometric phase of 
the particle flux mgf. 

To study diffusion without the external field, A(x, t) — 
0, Refs. [13] derived the path integral for the mgf by dis- 
cretizing the space into small intervals of length a <^ L, 
indexed by i. Then Poisson transition rates among the 
neighboring intervals are prescribed in a way that the 
continuous limit a — > recovers the diffusion equation. 
This reduces the path integral derivation to an already 
solved problem of stochastic transitions among a dis- 
crete set of states. To include the force A{x,t), we 
assume that it creates an asymmetry in the left and 
right transition rates. For example, (13 1 can be re- 



covered if the transition rates are such that during a 
short time 5t the average numbers of particles trans- 
fered left and right are (SQi^i^i) = Dp{xi)St/a and 
(SQi^i+i) — Dp{xi)5t/a + A{xi)p{xi)6t, respectively. 
Then repeating the same steps as in |13| and taking the 
continuous limit, we find the following path integral rep- 
resentation of the generating function: 

Z= ^ j Dp{x,t)Dxix,t)e-l'o^°'^*-l'o'^''^'^P+"^P^^^^ 



(14) 

where Q is the difference between the numbers of par- 
ticles passing through a; = in the right and the left 
directions during the period To- The Hamiltonian is 

iIte,).-M(.,,v| + iD||-D.(|)',(15) 



The dependence on the counting field xc in (15 1 is 
hidden in the boundary conditions on x [12 , which. 



for a periodic system with the spatial period L, are 
p{L) = p(0), and xW = x(0) - xc- Now, solving the 
saddle point equations and substituting the result back 
into the action in the path integral, we write the mgf 
in a famihar form Z{xc) = exp{5'geom(xc) + Sdixc)}, 
where Sg^omixc) = Jo° dt dx {iXciPd), and Sdixc) = 

In" '^^lo dxH{pci{xc),Xci{xc))- 

The analysis simplifies if we are interested only 
in mean currents, rather than in their fiuctuations. 
Then we consider xc ^ 1 and find the contribution 
to logZ that is linear in it. In fact, only 5*200111 
has this contribution in our case. To determine it, 
it is sufficient to find pci{x,t) to the zeroth order 
and Xci(a;,t) to the first order in xc- This results 
in poi(x,t) « [Qo/i?-(i)]e-^("^*)/'=«^, Xci{x,t) « 
[-Xc/i?+(i)]/o"e^("'^*)/'=«^da;', where i?±(t) 
jL^±v{.j)/ksT^^^ and Qo - jf p,,{x,t)\^^=odx is 
the number of particles per unit cell. This leads to 
Z{xc, Tq) = exp[ixcJTQ + O (xc)] ' where the terms in 
0{xc) *^an reveal the higher order cumulants, and the 
average current J — ~{i/TQ){d^^ log Z)^^^q is 



J 



2Tf 



Jo 



To pL 

dt / dx {dtv dxU — dxV dtu) , 
Jo 



(16) 



where we introduced u{x,t) 



R-{t) Jo 



j-g-y(.',t)/fe«T^^, 



a.ndv{x,t) = T^j„ 
erty u{L, t) = v{L, t) ~ \. The integrand in ( 16 1 is a pure 
curl of a vector A with components Ax = v{x, t)dxu{x, t) 
and At = v{x,t)dtu(x,t) defined in the two dimensional 
space-time. Thus the current can be expressed as 



J 



2To 



A dr, 



(17) 



where dr = {dx, dt), and c is the contour that encloses a 
space-time cell with boundaries a,t x — 0, L and t — 0,Tq. 

For a uniformly shifting potential V{x, t) — V{x — 
tL/Tp), R± are time-independent, and the integration 
in ^ leads to J = Qo/Tq - {Qo/To)L^ /{R+R_). The 
first term in this expression is the quantized contribution 
which is dominating in the limit of a large potential am- 
plitude. In [3], this quantization of the classical ratchet 
current was connected to the Chern number of the Bloch 
band related to the potential V{x). 

Pump current in the SIS epidemiological model. In a 
final calculation, we show how the stochastic path inte- 
gral allows derivation of pump-like effects in novel sce- 
narios; specifically where, unlike in our first example, 
the system cannot be factored into non-interacting iden- 
tical stochastic subsystems. We consider the standard 
Susceptible-Infected-Susceptible (SIS) mechanism of an 
infection outbreak, which is a good model for infiuenza. 
State of the art epidemiological modeling uses determin- 
istic dynamics [161, which tracks only fractions of popu- 
lations in various states during an outbreak progression. 
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However, it is understood that stochasticity may be es- 
sential. Thus here we discuss if stochasticity, and espe- 
cially effects due to slow variability of the infectivity and 
the recovery rates, can affect disease outbreaks. 

Let's denote infected individuals by / and their number 
by N. The disease spreads due to a permanent infection 
source and because it can be transmitted by the infected 
individuals. All infected people eventually recover. 

Thus the full kinetic scheme is 

1. ^ /; rate fci (permanent infection source); 

2. / — > 0; rate per infected individual fc2 (recovery); 

3. / ^ I + I; rate per infected individual (infection 
spread by contacts). 

Here ki are independent of N because we assume that 
outbreaks are small in comparison to the total population 
size (still 3> 1 is possible). This requires that ^2 > fcs, 
so that, if stochasticity is unimportant, the determin- 
istic steady state solution is A^st = ^1/(^2 — ^3), and 
the stationary flux into and out of the infected state is 
Jst = k2Nst — fcifc2/(fe— ^3)- This model is a birth-death 
process, and, with time-independent rates, it has been 
studied extensively. Here we are interested in estimat- 
ing (possibly substantial) effects of rate time-dependence. 
The Hamiltonian in the path integral for this model is 

H{x, A, t) = ki{t)e^^ + k2Ne(^+^^) + kz{t)Ne^^, (18) 

where x is the conjugated variable to A^ and xc counts 
the flux out of /. With A" ^ 1, we can use the saddle 
point analysis, which is exact since H is linear in A. 

Now consider a periodic time dependence of the rates 
ki, which may be due to the time-of-day or seasonal ef- 
fects. For simplicity, we assume that only fci and k^ vary, 
and the recovery rate ^2 remains constant. As before, 
the mgf has both the classical and the geometric contri- 
butions, i.e. Z = exp[5'ci -t- ^goom]- The classical one is 
the average of the stationary mgf over the period of the 
rates variation, Tq, while the geometric one is again an 
integral over the surface Sc inside the contour enclosed 
by fci(t) and k^{t): 



Sc\ + 5gcom I dtH[xci{t),N,i{t),t] 

^0 Jo 



T 



- 11 dkidk3Fki,h,(k), 
J JSc- 



k2{K^ - 2fc3exc - k) 



2 hi/^ 



H{xchN,i,t) 



ki{K^ - k) 
2fa ' 



(19) 
(20) 
(21) 



where now 



This corresponds to the mean flux J = J, 



k^, and k = \ K"^ —4/02^36 



Jo 



Tn //s 



'^3<--xc ■ 

pump r Jci, where 



geometric contribution, and the classical flux is Jd = 
1/To Jstdt. Notice, in particular, that Jpump esc KZ^ , 
and it can become very large near K- = 0. Fluctuations 



are easy to compute as well by differentiating (19 1. 

Conclusion. Based on the stochastic path integral 
technique, we built the theory of geometric fluxes in clas- 
sical stochastic kinetics, and we proposed a general ap- 
proach for identiflcation and calculation of pump-like cur- 
rents, including the familiar reversible ratchet, as well as 
new phenomena. In the adiabatic limit, the full count- 
ing statistics of pump fluxes is provided by the term that 
depends on the choice of the contour in the parameter 
space, but does not depend on the rate of the motion 
along this contour, and thus has a geometric origin. The 
stationary saddle point approximation of the path inte- 
gral is sufficient for calculations of this geometric contri- 
bution in the case of a large number of particles. This 
approach leads to the complete theory of reversible effects 
in noncquilibrium statistical physics. It will open doors 
to a study of such poorly understood systems as ratch- 
ets with interacting diffusing particles, or epidemiological 
models on complex social networks with time-dependent 
parameters. 
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